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ABSTRACT 

Time-resolved Ka spectra of magnetically-active interacting binaries have been used to create 3D 
Doppler tomograms by means of the Radioastronomical Approach. This is the first 3D reconstruction 
of j3 Per, with RS Vul for comparison. These 3D tomograms have revealed evidence of the mass 
transfer process (gas stream, circumprimary emission, localized region, absorption zone), as well as 
loop prominences and coronal mass ejections (CMEs) in j3 Per and RS Vul that could not be discovered 
from 2D tomograms alone. The gas stream in both binaries may have been deflected beyond the central 
plane by the mass losers magnetic field. The stream was more elongated along the predicted trajectory 
in RS Vul than in j3 Per, but not as pronounced as in U CrB (stream-state). The loop prominence 
reached maximum V z velocities of ±155 km s _1 in RS Vul compared to ±120 km s _1 in /? Per, while 
the CME reached a maximum V z velocity of ±150 km s _1 in RS Vul and ±100 km s _1 in /? Per. 
The 3D tomograms show that the gas flows are not symmetric relative to the central plane and are 
not confined to that plane; a result confirmed by recent 15GHz VLBI radio images of f3 Per. Both 
the 3D Ha tomography and the VLBI radio images support an earlier prediction of the superhump 
phenomenon in j3 Per: that the gas between the stars is threaded with a magnetic field even though 
the hot B8V mass gaining star is not known to have a magnetic field. 

Subject headings: techniques: image processing - accretion, accretion disks - stars: binaries: close - 
(stars:) binaries: eclipsing - (stars:) circumstellar matter - stars: imaging - stars: 
individual (/? Persei, Algol, RS Vulpeculae) 



1. INTRODUCTION 

The evolution of interacting binaries occurs as a con- 
sequence of mass transfer between the stars, and this 
process has been studied using spectroscopic and pho- 
tometric analyses, theoretical simulations, and imaging 
techniques. Specifically, the image reconstruction tech- 
nique of Doppler tomography has provided images show- 
ing how the gas moves (or flows) between and around 
the stars. Moreover, tomography of time-resolved spec- 
tra has been used to determine observational measure- 
ments of the mass transfer rate in these binaries and 
to confirm the rates derived wi th the aid of the shell- 
spec spectrum sy nthesis code (|Budaj fc Richard's! 120041 : 
Mille r et al.l 12007). This derived rate is equal to the ini- 
tial mass loss rate at the onset of Roche Lobe overflow 
since the donor star will transfer gas to its co mpanion un - 
til the initial mass ratio has been reversed (|Iben| [l991). 
Hence, tomographic images of interacting binaries at var- 
ious evolutionary stages yield important clues about the 
life cycles of interacting binaries. The early to late evolu- 
tionary stages of interacting binaries have been studied 
through spectroscopic analysis, calculation of synthetic 
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spectra, and w ith both 2D and 3D tomography (e.g., see 
iRichardsl [20121 ). 

Two-dimensional tomography has provided images of: 
(1) accretion disks around the mass gaining star in the 
long-period Algols, CVs, nova-likes, X-ray binaries, and 
gamma-ray binaries; (2) gas streams and gas flowing 
along magnetic field lines in the magnetic CVs; (3) a com- 
bination of Keplerian disks, gas streams flowing along 
the predicted gravitational path, shock regions where the 
stream and disk interact, chromospheres, and other mag- 
netic structures in Algol binaries; (4) regions where the 
gas slows down after circling the mass gainer in direct- 
impact systems (e.g., Algol, j3 Per); and (5) asymmetric 
accretion disks that maintain their asymmetry on long 
timescales (e.g., TT Hya) . For a summary of the 2D 
tomography, see Richards ([20041 ). 

Three-dimensional Doppler tomography was intro- 
duced usin g the Rad i oastron omica l Approach (RA) de- 
veloped bv lAgafonovl (2004a. b) and Agafo nov fc Sharoval 
(2005a, b). A summary of the important factors 
in the realization of 3D tomography is given in 
lAgafonov fc Sharoval ([2013D . To date, the 3D RA tech- 
nique has been applied to only four systems: three Algol- 
type eclipsing binaries (U CrB, RS Vul, and now f3 Per) 
and the X-ray binary Cyg X-l. The first 3D tomo- 
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grams f or the entire class of interacting binarie s were cre- 
ated by Agaf onov. Richards fc Sharoval (j2QQ6h using the 
U CrB alternating Algol system, which displays disklike 
emission in the tomogram at some epochs and stream- 
like emission at other epochs. The 3D tomograms of 
th e disk and stream states of U CrB were described 
by Agafonov, Sharova & Richards (2009) and compared 
to RS Vul by iRichards. Sharova fc Agafonovl (|201Q[ ). A 
prominent gas stream was also discovered in the 2D and 
3D images of Cyg X-l (jSharova et al.ll2012l ). 

The 3D images illustrate the distribution of gas flows 
in the orbital plane and beyond that plane, and confirm 
the range and complexity of the emission sources in these 
binaries: from accretion disks and gas streams, to shock 
regions and magnetic structures. Additional information 
provided by 3D tomography includes evidence that (6) 
prominent and extensive gas flows exist beyond the cen- 
tral plane of the binary in the z direction (e.g., U CrB, 
RS Vul); (7) the accretion disk may precess or may be 
tilted relative to the central plane (e.g., U CrB); and (8) 
loop prominences and coronal mass ejections associated 
with the magnetic field of the donor star also contribute 
to the gas flows (e.g. , RS V ul). These results were sum- 
marized by IRichards] (|2012l ) . 

In this paper, we continue our exploration of 3D tomog- 
raphy by comparing two magnetically-active systems: j3 
Per (Algol) and RS Vul. Specifically, we have 1) identi- 
fied the features in the 3D images and related them to 
those seen in the 2D tomograms; 2) compared the 3D 
tomogram of /? Per with those of RS Vul derived earlier; 
and 3) provided a physical interpretation of these fea- 
tures. The system parameters are described in Section 
2, the 2D tomograms are reviewed in Section 3, the 3D 
tomograms are described in Section 4, a model for /? Per 
and RS Vul based on the 3D images is given in Section 
5, and the conclusions are stated in Section 6. 

2. SYSTEM PROPERTIES FOR PER AND RS VUL 

/3 Per (Algol) is the prototype of the Algol-type bi- 
naries. It is the brightest and nearest eclipsing binary 
(d=29 pc), and one of the most extensively studied ob- 
jects in the sky at all wavelengths. It contains a partially 
eclipsing binary , with a tertiary in a 1.86 ye ar orbit about 
the binary (see Richa rds fc Albrightl IT999). The system 
consists of a B8V primary, a K2IV secon dary, and an 
F1V tertiary ([MorgaiJ 119351 : IRichards! fl993h . 6 Per is a 
short-period Algol (P orb = 2.8673 days; Iffill et~a71 Il971h 
with a mass ratio, g=0.22 ±0.03 (Richards 1993) an d 
an orbital inclination of 81. 4±0. 2° (Richards et al. 1988). 
Other system properties include: M p = 3.7±0.3 Mq, M s 
= 0.81 ± 0.05 Mm, Rn = 2.90 ±0.0.04 R A = 3.5 ±0.1 
R Q (jTomkin fc Lambertlfl978l; [Ric hards 1993), sy stemi c 
velocity V Q = 3.8 km s" 1 ([Hill et all 11975 lEatonlll975f ). 
and velocity semi- amplitude K v = 44.0 ±0.4 km s -1 
([Hill et al.lll97ll : lTomkin fc L amber til 19 78[ ). iBaron et al.1 
(2012) used observations from the CHARA interferome- 
ter to determine new masses, radii and orbital properties 
of /? Per that are consistent with the parameters derived 
earlier. 

The properties of Per are compared with those of RS 
Vul in Table 1; see IRichards. Sharova fc Agafonovl ([20101 ) 
for the RS Vul references. One notable difference is that 
RS Vul contains a Gl secondary compared to the K2 
mass loser in f3 Per. Figure 1 also shows the similarities 



TABLE 1 
Properties of (3 Per and RS Vul 



Property (3 Per RS Vul 



Spectral types 


B8V+K2IV+F1V 


B5V+G1III 


P orb (days) 


2.867315 


4.4776635 


mass ratio, q 


0.22 ±0.03 


0.27 


in 


81.4±0.2° 


78.7±0.2° 


M p (Mq) 


3.7±0.3 


6.59±0.15 


M s (Mq) 


0.81 ±0.05 


1.76 ±0.05 


R P (Rq) 


2.90 ±0.0.04 


4.71 ±0.48 


Rs (Rq) 


3.5 ±0.1 


5.84 ±0.23 


V (km s- 1 ) 


3.8 


-20.1 


K p (km s- 1 ) 


44.0 ±0.4 


54.0 ±1.0 



between the Roche geometries of the two binaries. Both 
binaries have similar mass ratios and the R p /R s ratio is 
~ 0.8 for both systems even though the stars in RS Vul 
are larger than those in Algol, hence these binaries are 
expected to have similar accretion structures. 

Magnetic activi ty associated with Algol se condaries is 
well known (e.g., Richar d^ fc Albrightl 1 19931 ). and both 
f3 Per and RS Vul contain late-type magnetically-active 
secondaries. /3 Per is a strong radio and X-ray source 
(Wade fe Hiellmindll97l iWhite et alJUoSfl IStern et al.l 
1992. Il995l ). and its cool secondary star is the dominant 
source of the emission at those wa velengths, character- 
istic of gyrosynchrotron radiation (jvan den Oord et al.l 
119891 : iFavata et al"]l2QQQD . Radio emission from RS Vul 
was detected with a flux density of 0.26 mJy at 5 GHz 
(jUmana et gil.l Il998h compared to flux densities as high 
as 3 mJy at 2.3 GHz a nd 1.2 mJy at 8.3 GHz for f3 
Per (Richards et al. 2003). RS Vul also has an X-ray lu- 
minosity of 2.0 x 10 30 erg s -1 ([White fc Marshallll 19831 ). 
that is about half of Algol's luminosity. So, the magnetic 
field of Algol may be more the active than that of RS Vul. 
Hence a direct comparison between the 3D tomography 
and the new radio image should be interesting. 

3. 2D DOPPLER IMAGES OF (3 PER AND RS VUL 

Two-dimensional Doppler tomograms of f3 Per were 
created from one hundred and thirty-five time-resolved 
Ha (6562.8 A) spectra of /3 Per, with a dispersion of 
0.093 A/pixel (4.3 km s _1 /pixel) or 6.9 Amm" 1 , col- 
lected with the solar-stellar spectrograph on the 1.5m 
McMath-Pierce Telescope at the Na tional Solar Observa- 
tory from 1992 October 6-21 (see Richar ds fc Alb right 
1999). The spectra were obtained at closely spaced po- 
sitions around the entire orbit of the binary, and all of 
the spectra used with the tomography procedure were 
obtained at phases outside of primary eclipse to exclude 
any contributions from the Rossiter-McLaughlin effect. 
Since the observed spectra are dominated by the spec- 
trum of the B8 mass gaining star, difference spectra were 
calculated by subtracting theoretic al photospheric spec- 
tra fr om the observed spectra (see IRichards fc Albrightl 
119991 for details). These difference spectra were used to 
calculate the tomograms of j3 Per. 

Figure 1 shows the 2D Doppler tomogram based on the 
1992 spec tra of /3 Per alon g with the 2D tomogram of RS 
Vul (from lRrchardsll200lh based on Ha spectra collected 
in two observing runs over 12 days from 1993 Apr 30 - 
May 31. These 2D images are later compared with the 
3D Doppler tomograms based on the same data. In the 
tomograms, the solid trajectory is the gravitational free- 
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Fig. 1. — Two-dimensional Ha Doppler tomograms and contour maps of /3 Per - 1992 (upper middle) and RS Vul - 1993 (lower middle), 
with Cartesian models of the two systems (l eft) and interpol ated trailed spectrograms (right). In the tomograms, the brightest sources are 
in deep red and the faintest in blue. Credit: Richards (2012). 



fall path of the gas stream; and the circles along this tra- 
jectory are marked at intervals of a tenth of the distance 
from the Li point to the distance of closest approach to 
the mass gainer. The largest solid circle and the smaller 
dashed circle mark the inner and outer edge of a Ke- 
plerian disk, respectively; and the plus sign marks the 
center of mass of the binary, which is close to the stellar 
photosphere in the cases of j3 Per and RS Vul. The veloc- 
ity images show the center of mass at (V x ,V y ) = 0, while 
the secondary star is placed at the velocity correspond- 
ing to its velocity semi-amplitude. In the tomograms, 
the brightest sources are in deep red and the faintest in 
blue. 

The first 2D tomogra m of Per, based o n 1976/77 
spectra, was created by Richa rds et al.l (jl996[ ). That to- 
mogram is similar to those derived from spectra collected 
in 1992 and again in 1994, more than 16 years later. The 
1976/77 trailed spectra of f3 Per displayed both single- 
and double-peaked emission whose orbital variations fol- 
lowed S-wave patterns. The sources of emission identi- 
fied in the 2D tomogram were, from strongest to weakest: 
(1) a gas stream: the source of the single-peaked emis- 
sion with velocities near the L\ point; (2) a localized 
region: concentrated in the region directly between the 
stars and along the line of centers; (3) a sub-Keplerian 
crescent-shaped disk-like structure: the source of the 
double-peaked emission; and (4) an emission source on 
the velocity of the secondary star provided evidence of 
magnetic activity associated with that star. An indepen- 
dent S-wave analysis performed by Ric hards et al.l {l996) 
found that the source of the single-peaked emission was 
coincident with the part of the gas stream near the L\ 
point, while the source of the double-peaked emission 



was coincident with the concentrated Localized Region 
found in the 2D image. 

4. THE 3D TOMOGRAPHY IMAGE RECONSTRUCTION 

The extension from 2D to 3D tomography has 
provided unprecedented views of the gas flows be- 
yond the central plane in which the stars orbit. In 
this work, we u sed th e Radioast ronomical Approach 
of lAgafonovl (|2004al lbh and lAgafonov fc Sharoval 
(2005a. b) that was dev eloped further by application 

to four binaries fsee lAgafonov. Richards fc Sharoval 

2006L lAgafonov. Sharova fc Richardsl I2009L and 
Richards. Sharova fc Agafonovl l2010f ). A summary 
of the important fac tors in the realization of 3D 
tomography is given by Agafo nov fc Sharoval (j2013f ). 

The basic assumpt ions and constraints of the 3D tech- 
nique include (jRichards. Sharova fc Agafonovl l2010): (1) 
the spectra are assumed to be broadened primarily by 
Doppler motions; (2) spectra with high wavelength res- 
olution and good resolution in orbital phase (or projec- 
tions) are essential; (3) spectra dominated by emission 
lines are primarily used; (4) the assumption that the 
gas is optically thin is sufficient to create a good first- 
order view of the binary; (5) the orbital inclination influ- 
ences the velocity resolution in the V z direction relative 
to those in the V x and V y directions, hence an adjust- 
ment is needed to compensate for the stretching effect 
of the Summa rized Point Spread Function (SPSF) i n the 
V z direction (Ag afonov. Richards fc Sharoval 12006): (6) 
the transformation of the emission intensity from veloc- 
ity space to coordinate Cartesian space is very compli- 
cated and it has not yet been solved even for 2D images. 
The RA technique creates 3D velocity maps of the gas 
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Fig. 2. — Comparison between the original data (left frame) and the spectra computed from the reconstructed 3D Doppler map of /3 Per 
(middle frame) in terms of the radial velocity, V r versus orbital phase. The right frame displays the orbital phase variation of the relative 
chi-square statistic, x 2 /Xc> where Xc * s the critical value corresponding to the 99% confidence level. The agreement between the observed 
and computed spectra is very good. 



intensity and it is an efficient method of processing the 
image in the case when only a limited number of spectra 
are available. 

4.1. The 3D Tomogram of (3 Per 

The 3D tomogram of /? Per was constructed from the 
same Ka spectra from 1992 October that were used 
to pr oduce the 2D Doppler tomogram (see Ric hards! 
I2QQU ). The 3D image was calculated with dimensions 
(Vx,V y ,V z ) for values ranging from -800 to +800 km s _1 , 
with intensities normalized to permit a comparison be- 
tween slices. The velocity resolution of the reconstructed 
3D image depends on the number of projections (i.e., 
number of spectra) and their distribution in orbital 
phase; while the orbital inclination influences the ratio of 
the resolutions in the V x , V y: and V z directions. The 4.3 
km s _1 dispersion is high and certainly good enough to 
permit an optimal resolution of the 3D image since the 
dispersion plays a smaller role in setting the resolution of 
the image compared to the constraints set by the orbital 
inclination. 

Figure 2 shows the comparison between the observed 
spectra of j3 Per and those computed from the recon- 
structed 3D tomograms. This figure was used to exam- 
ine the quality of the reconstructed 3D image produced 
by the RA method and it shows the trailed spectrograms 
in which the radial velocity, V r is plotted versus orbital 
phase. The observed spectra are displayed in the left 
frames and the spectra computed from the 3D tomo- 
grams are displayed in the middle frames. Two S-wave 
patterns are noticeable in the observed and computed 
spectra corresponding to a broad emission component 
and a narrow emission component. The overall agree- 
ment between the observed and computed spectra is very 
good. The quality of the fit is shown in the right frame of 
Figure 2, which displays the relative chi-square statistic, 
X 2 /Xc> versus orbital phase. Here, x 2 is normalized to 
the critical value of x 2 , which corresponds to the largest 
acceptable value of x 2 at the 99% confidence level and 
demonstrates the quality of the calculated values of x 2 - 
Figure 2 shows that x 2 was l ess than 30% of the critical 
value at most orbital phases and it was only as high as 
50% of the critical value over a narrow phase range from 



0.15 to 0.20, perhaps because the 3D code assumes that 
the gas is optically thin when it might be optically thick 
at these phases. Since x 2 was much lower than the crit- 
ical value the 3D reconstruction provides a good match 
to the data. 

The 3D tomogram of j3 Per is displayed in Figure 3 as 
fifteen slices in the (V x ,V y ) plane from V z = —420 km s _1 
to +420 km s _1 , at intervals of 60 km s _1 . Two addi- 
tional slices in the (V y ,V z ) direction are shown in the 
bottom frames. Figure 3 shows that there are usually 
relatively small changes in the distribution of gas veloc- 
ities between slices for a V z grid spacing of 60 km s -1 . 
The most interesting representative slices in the 3D to- 
mograms of f3 Per and RS Vul are compared in Figure 4; 
for comparison, the full s et of slices in the 3D tomogram 
of RS Vul can be found in lRichards. Sharova fc Agafonov! 
([2010h . 

Since the orbital inclination of j3 Per is high (i = 81.4°), 
the velocity resolutions in the V x and V y directions will 
be better than the resolution in the V z direction. The 
resolution in the V z direction is the same as that in the 
other directions only for an inclination of 45°, and de- 
grades as the inclination increases toward 90°. Conse- 
quently, the 3D tomogram of j3 Per was restored with a 
resolution of 40 km s _1 in the V x and V y directions and 
180 km s _1 in the V z direction, corresponding to a SPSF 
half power beam width (HPBW) of 40 x 40 x 180 in units 
of km s _1 . In order to adjust for the artificial stretch- 
ing effect of the SPSF, the V z velocities were scaled by a 
factor of 0.22 (= 40/180), corresponding to the ratio of 
the velocity resolutions in the V z direction compared to 
the other directions. This uniform adjustment is useful 
in illustrating the locations of the emission sources in the 
3D image. The adjusted V z values are given in the last 
column of Table 2. 

Figure 5 illustrates how the SPSF affects the resolution 
of the image. Its pattern on the HPBW in the (V x y y ) 
plane (see left frame) leads to the equal resolution along 
the V x and V y axes. However, the image is also stretched 
approximately four times in the V z direction relative to 
the other directions (see the (V y ,V z ) plane in the right 
frame). An example of the deconvolution to adjust for 
this artificial distortion is illustrated in Figure 6. This 
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Fig. 3. — Visualization of the (V x ,V y ) 2D slices in the 3D Doppler Tomogram of (3 Per (1992) displayed symmetrically from V z = -420 
km s _1 to V z = +420 km s _1 , in steps of 60 km s _1 in the V z direction. Two additional images in the (V y ,V z ) direction are shown (bottom 
frames). 
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Fig. 4. — A comparison of symmetric (V x ,V y ) slices along the V z direction in the 3D Doppler tomograms of ft Per and RS Vul. The blue 
dashed lines illustrate the trajectory of the gas that flowed around the mass gainer only to be slowed by impact with the incoming gas 
stream and the Localized Region (shown for both binaries). 



TABLE 2 

Characteristics and locations of prominent emission features in the 3D tomogram of (3 Per 





Location: central velocity or velocity range 


(km s- 1 ) 


No. Emission Feature 

i rAi • • • 


V X Vy V Z 


V z (adjusted) 



2 Emission on donor star (active regions) 
Emission near LI 

3 Gas stream flow 

4 Star-stream impact region 

5 Locus of the accretion disk 

6 Localized Region (LR) - Part 1 

7 Localized Region - Part 2 

8 High velocity flow near donor 

9 Absorption zone 



(-70 to +70) 180 (80 to 280) 

+110 

-520 to a to 300 a 
No source found 
No source found 

(-50 to +50) +50 (0 to 100) 

70 (20 to 120) +50 (-20 to 140) 

110 (50 to 170) 220 (150 to 300) 

(-60 to +60) -60 (0 to -120) 



(-500 to +500) 
-0 (-500 to +500) 
-300 (-480 to +120) 



200 (0 to 400) 
60 (-180 to 240) 
340 (240 to 440) 
240 (80 to 460) 



(-120 to +120) 
-0 (-120 to +120) 
-70 (-115 to +30) 



50 (0 to 100) 
15 (-45 to 60) 
80 (60 to 100) 
60 (20 to 115) 



a corresponds to the predicted ballistic trajectory of the gas stream. 
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Fig. 5. — Slices of the SPSF calculated for all the 3D viewing 
angles (inclination and orbital phases) of the (5 Per binary. The 
SPSF affects the resolution of the image in a symmetric way in the 
(y x ,V y ) plane (left frame) while artificially distorting the image 
along the V z direction in the (V y ,V z ) plane (right frame). This 
effect is strongly influenced by the inclination of the binary to the 
plane of the sky. 
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Fig. 6. — Cross sections of the (3 Per and RS Vul 3D Doppler 
tomograms in the (V y ,V z ) plane for V x =0 km s -1 . The top frames 
show the original images, and the bottom frames show the cen- 
tral parts of the same slices after an approximate deconvolution 
to adjust for the artificial stretching effect of the SPSF in the V z 
direction. The numbered features correspond to those in Figs. 3 
and 4. 

figure shows cross sections of the 3D tomograms in the 
(V y: V z ) plane for V x =0 km s _1 , before and after the ad- 
justment to the velocities in the V z direction. A uniform 
scaling of the images by the factor of 0.22 provides a 
rough first approximation to the complex deconvolution 
that is needed to rectify the images. 

4.2. Description and Interpretation of Emission 
Features 

The dynamical properties of the mass transfer process 
in j3 Per were expected to be similar to those of RS Vul 
based on their Roche lobe geometries shown in Figure 1. 
The 3D tomograms shown in Figures 3 and 4 confirm that 
the 3D reconstruction has reproduced the 2D tomograms 
for the case when V z = km s _1 . Moreover Figure 4 
shows some general similarities between the 3D images 
of these two binaries. 

Several main features and other weaker features were 
detected in the tomograms shown in Figures 3 and 4. The 
numbering system for these features is given in Table 2 



and is consistent with earlier papers in this series. The 
V z velocities, V z (adj), quoted below are the values after 
the adjustment for the stretching effect in this direction. 

The strongest sources in f3 Per are listed below. 
Feature 1: An emission source associated with the pri- 
mary star, the circumprimary emission, was detected 
with high negative V z (adj) velocities from -50 to -150 
km s _1 in /? Per compared to to +90 km s _1 in RS Vul. 
This source provides evidence that the mass gainer was 
spun- up to super-synchronous values by the impact of the 
gas stream onto the stellar surface. Moreover, the direc- 
tion of the flow in the V z direction in j3 Per was opposite 
to that in RS Vul. In both binaries, the projected rota- 
tional velocity, vsini is comparable to the synchronous 
rotational velocity of the mass gainer (50 - 60 km s _1 ), 
and both are lower than the maximum V z (adj) veloci- 
ties associated wi th this feature (e.g., vsmi = 53 ± 3 
km s -1 for f3 Per; iRucinskil [1971 iTomkin fc fanlll985f ). 
Since the V z (adj) velocities associated with Feature 1 are 
higher than the star's rotational velocity, it is reasonable 
to assume that the star may have been spun up, pre- 
sumably by the impact of the gas stream onto the stellar 
photosphere. Similar super-synchronous rotation veloci- 
ties which have been measured directly from line profile 
analysis of direc t-impact systems including f3 Per (e.g., 
lMallamalll978allbh . Velocities of ~500 km s" 1 at the 
impact site were predicted by iLubow fc Shul (|1975f ) us- 
ing semi-analytical ballistic calculations of direct-impact 
systems; such high velocities could produce supersonic 
turbulence and increase the angu lar momentum of the 
material around the mass gainer (|Richardslll993l ). 
Feature 2: An emission source associated with active 
magnetic regions on the donor-star was detected with 
a wide range of V z (adj) velocities from -120 to +120 
km s _1 , equally distributed about the central plane. 
In the (V x , V y ) slices for both binaries, this feature has 
nearly the same velocity as the donor star. However, in 
RS Vul, this emission displayed a wider V z range from - 
150 to +155 km s _1 , suggesting higher magnetic activity 
than found for /? Per. This behavior over a 305 km s _1 
range in V z (adj) led to the discovery of a loop promi- 
nence in RS Vul along which gas is rising and falling since 
the footpoint in the (V x , Vu) plane stayed the same while 
the V z velocity changed ([Richards, Sharova fc Agafonovl 
120101 ). The V z (adj) range of 240 km s _1 in /? Per suggests 
that yet another prominence has been found in Algol- 
type binaries. 

Feature 3: An emission source associated with the gas 
stream along the predicted gravitational trajectory in the 
(V X: V y ) plane was detected at V z (adj) velocities from - 
115 to +30 km s _1 . This source reached the same high 
velocities found in RS Vul, but it was not as prominent 
(or bright) in j3 Per as in RS Vul. In both binaries, this 
source reached peak extension at high negative V z ve- 
locities, and not in the central plane (V z = km s _1 ), 
suggesting that the gas stream may have been deflected 
away from the central (V X: V y ) plane by the magnetic ef- 
fects of the mass losing stars in both binaries. 
Feature 4: The star-stream impact region was not de- 
tected in either binary even though the gas stream flow 
reached the velocity associated with the surface of the 
mass gaining star in both cases. 
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Feature 5: Very little emission was found within the 
locus of a Keplerian accretion disk in Algol, and only a 
small amount of gas was detected in the case of RS Vul. 
Features 6,7: The Localized Region (LR) between the 
stars was detected in both binaries over a range of V z 
velocities. This source represents the part of the gas flow 
that circles the star and interacts with the incoming gas 
stream. The blue dashed lines in Figure 4 illustrate the 
trajectory of the gas that flowed around the mass gainer 
only to be slowed by impact with the incoming gas stream 
and the Localized Region. 

Feature 8: The high velocity flow associated with the 
donor star reached V z (adj) velocities of +100 km s _1 in 
j3 Per, compared to +150 km s _1 in RS Vul. Relative to 
the general accretion flow, these large positive velocities 
suggest that gas was ejected from the donor star in a 
manner similar to a coronal mass ejection (CME), which 
can achieve velocities from hundreds of km s _1 to over 
1000 km s _1 in extreme cases on the Sun. This result 
suggests once again that the cool donor star was more 
active in RS Vul than in /? Per. 

Feature 9: An absorption source associated with the 
velocity of the mass gaining star in the (V x ,V y ) plane 
was detected at positive V z (adj) velocities, from +20 to 
+115 km s _1 . In RS Vul, this feature had higher veloc- 
ities in the (V x ,V y ) plane, over a wider range of veloci- 
ties in that plane, and for V z (adj) velocities from -35 to 
+155 km s _1 . This source provides additional evidence 
that the mass gainer was spun-up by the impact of the 
gas stream onto the stellar surface and heated to temper- 
atures beyond the optical regime. A similar effect was 
found for U Sge when the Ka tomogram was com pared 
with the uv tomogram ([Kempner fc Richar ds 1999). The 
uv tomogram showed that the hotter gas was roughly in 
the same part of the tomogram as the absorption region 
near the accreting star. 

5. A MODEL FOR THE MAGNETICALLY ACTIVE 
BINARIES OF /3 PER AND RS VUL 

In summary, the general features found in 2D slices 
of the 3D tomograms of the interacting binaries U CrB, 
RS Vul, f3 Per, and Cyg X-l include: (1) circumprimary 
emission, (2) emission from active magnetic regions asso- 
ciated with the donor star, (3) flow along the predicted 
trajectory of the gas stream, (4) the star-stream impact 
site, (5) the predicted locus of a Keplerian accretion disk, 
(6) the localized region between stars, (7) a second local- 
ized region, (8) high velocity flow moving away from the 
donor star, and (9) the absorption zone near the mass 
gaining star. The characteristics of most of these sources 
found in /3 Per were eithe r identical or similar to those 
found for RS Vul (see iRichards. Sharova fc Agafonov! 
l2010f ). However, the star-stream interaction (4) was ab- 
sent in both binaries, and the Keplerian disk (5) was 
absent in j3 Per and weak in RS Vul. The absorption 
zone (9) in f3 Per was found with low velocities near the 
mass gainer, but with a higher range of velocities in RS 
Vul. Figure 7 shows the 3D representations of the ve- 
locity distribution of emission sources relative to the V z 
axis for both binaries. 

The emission sources associated with the cool magnet- 
ically active star (Features 2 and 8) are of particular 
interest in this work. Both j3 Per and RS Vul displayed 



evidence of loop prominences (2) and coronal mass ejec- 
tions (8). Moreover, these flows are not symmetric rel- 
ative to the central plane. Since the tomogram reveals 
the Doppler motion of the gas, then the detection of a 
smooth transition from positive V z to negative V z veloc- 
ities relative to a fixed footpoint in the central (V x ,V y ) 
plane is evidence of a smooth motion toward and away 
from the observer in the V z direction beyond the central 
plane (see Table 2). This behavior is characteristic of a 
loop prominence similar to those detected on the Sun. 
Similarly, a smooth transition from low to high positive 
V z velocities relative to a fixed footpoint in the central 
(V X: V y ) plane suggests motion away from the observer at 
high velocities. This behavior is reminiscent of the mass 
loss associated with flares known as CMEs. For both 
phemenona, RS Vul achieved higher V z velocities than j3 
Per. Moreover, 2D Doppler tomograms based on spec- 
tra of j3 Per collected from 1976 to 1994, showed that 
the image of /3 Per remaine d almost the same ov er that 
18-year period of time (e.g.. IRichards et aT1ll996[ ). 

Both of these phenomena have now been detected in 
/3 Per and RS Vul with the help of 3D Doppler tomog- 
raphy. This behavior was to be expected based on the 
radio and X-ray evidence that both systems should be 
magnetically active (See Section 2). A 5 . 6-year continu- 
ous radio flare survey by IRichards et al.l ([2003) revealed 
that radio flares occur regularly every 48.9 ± 1.7 days on 
the cool mass- losing star in j3 Per at both 2.3 GHz and 
8.3 GHz. This cycle corresponds to a major flare every 
17 orbital periods or roughly every 1.6 months. Hence it 
should not be surprising that a CME was detected in a 
tomogram collected over a period of 16 days in October 
1992. RS Vul was too faint to be included in that survey, 
so its flaring period is still unknown. 

The continuous radio survey of f3 Per at 2.3 GHz 
and 8.3 GHz also yielded evidence of superhump-like 
behavior, suggesting that the gas stream from the 
cool star is threaded by th at star's magnetic field 
([Retter. Richards fc W u 2005). This was the first detec- 
tion of the superhump phenomenon in the radio regime 
and the first observation of superhumps in Algol sys- 
tems. Moreover, this result suggests that the gas in 
the vicinity of the mass gaining star (flowing along the 
path of a sub-Keplerian annulus) in j3 Per is threaded 
with a magnetic field even though this B8V star is not 
known to have a magnetic field. Instead, the magnetic 
field in the annulus originates from amplification of the 
seed field in the magnetized material transferred from 
the mass losing star. Since the ram pressure of the ac- 
cretion flow is larger than the magnetic stress at the La- 
grangian p oint, the ionized accreting; ma terial drags the 
field along ([Retter. Richards fc Wull2005l ). The field lines 
are wound around the mass gaining star through the gas 
stream; the field is stretched and amplified, and at some 
stage reconnection can occur, causing flaring and accel- 
eration of the particles. The exception in the case of 
/3 Per is that only one star in the system has a strong 
magnetic field. In spite of the field reconnection in the 
gas around the mass gainer, the lack of radio emission 
near the B8V star could be explained by the Razin effect 
and/or plasma low- frequency cu toff when the electron 
density is high (e.g.. iDulk 111985ft . 

A consequence of this superhump effect is that the cir- 
cumprimary gas is expected to precess. This is interest- 
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Fig. 7. — 3D representations of the velocity distribution of emission sources relative to the V z axis for RS Vul and f$ Per. The top frames 
illustrate a tilted view of the velocity distribution of the gas beyond the central velocity plane of each binary. The other 3D views are 
shown in order of orbital phase starting with <f> = 0.0 (second row), <j) = 0.25 (third row), <f> = 0.5 (fourth row), and <f> = 0.75 (last row). 
The emission sources in the image have been scaled according to intensity, with the brightest sources shown as large red points and the 
fainter sources as small green points. 



mg since lAgafonov. Sharova fc Richards! (|2Q09l ) used 3D 
Doppler tomography to predict that the Keplerian accre- 
tion disk in another Algol- type binary, U CrB, does pre- 
cess. So, it is satisfying to discover that the super hump 
phenomenon is not restricted to the classical Keplerian 
accretion disks found in compact binaries. 

The superhump res ult predicted by 

Retter, Richards & Wu (2005) was confirmed five 
years later by iPeterson et al.l (|201Q[ ) through a global 
very long baseline radio interferometer array (VLBI) 
study which produced a resolved 15 GHz radio image 
of j3 Per based on observations collected in 2008. This 
radio image displayed evidence of gyrosynchrotron 
filled-loops associated with the cool donor star, a 
magnetically-threaded gas stream between the stars, 
and asymmetric magnetic structures that were not 
confined to the orbital plane; however, no radio emission 
was directly associated with the mass gaining star. 
Figure 8 shows the comparison between the VLBI 
results and those derived from 3D Doppler tomography 
of j3 Per. The resemblance suggests that tomography is 
truly a powerful tool that provides information about 
the gas flows that is consistent with the radio image of 
the binary, at least in the case of j3 Per. Our 3D results 



confirm the IPeterson et al.l (|2010f ) results, and also 
provide additional information about the interactions 
between the two stars beyond the orbital plane. 

6. CONCLUSIONS 

We have used time-resolved Ha spectra of two 
magnetically-active interacting binaries to extract infor- 
mation about the three-dimensional velocity distribution 
of the gas flows in these systems. This is the first 3D re- 
construction of f3 Per using the Radioastronomical Ap- 
proach. Emission and absorption features detected in 
2D tomograms have been confirmed in the 3D Doppler 
tomograms of j3 Per and RS Vul. The 3D images have 
revealed evidence of the mass transfer process, including 
the gas stream, circumprimary emission probably caused 
by the impact of the gas stream onto the surface of the 
hot mass gaining star, an absorption zone where gas has 
been heated to temperatures beyond the Ha regime, and 
a localized region where the gas that circled the mass 
gainer has been slowed to sub-Keplerian velocities by im- 
pact with the incoming gas stream. The 3D images also 
display evidence of magnetic activity associated with the 
cool star, including evidence of loop prominences and 
coronal mass ejections that could not be discovered from 
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15 GHz VLBI radio images (Courtesy: Peterson et al. 2010, Nature, 463 , 207) 
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Fig. 8. — Comparison between the 3D tomog raphy Ha i mag es of (3 Per (top frames) with the 15GHz VLBI radio images of the binary 
at two phases (lower frames) derived from [Peterson et al. ( 2010). The top images are velocity images while the lower images are position 
images. Magnetic activity can be detected both in the radio and Ha spectra and these images display consistent information about the 
binary. 



the 2D tomograms alone. 

The extension from 2D to 3D images has provided 
a view of the gas flowing beyond the central plane de- 
fined by the orbital motions of the stars. The 3D im- 
ages of f3 Per and RS Vul have once again revealed that 
the gas flowing beyond the central plane is substantial 
since the gas achieved V z (adj) velocities as high as ±155 
km s _1 . Specifically, the loop prominence reached max- 
imum V z (adj) velocities of -150 to +155 km s _1 in RS 
Vul compared to -120 to +120 km s _1 in /3 Per, and 
the CME reached a maximum V z (adj) velocity of +150 
km s _1 in RS Vul and +100 km s _1 in /3 Per. Hence, 
the gas velocities associated with magnetic activity on 
the cool star in RS Vul were slightly higher than those 
achieved by j3 Per. Moreover, the (V x ,V y ) slices showed 
that gas stream was more elongated in RS Vul than in 
/3 Per, but not as impressive as the distinctive extended 
flow along the predicted path seen in the 3D image of U 



CrB. In both systems, the stream path may have been 
deflected beyond the central plane by the magnetic field 
of the cool star. 

Both the 3D Ha tomography and the VLBI radio im- 
ages support an earlier prediction of the superhump phe- 
nomenon in j3 Per: that the gas between the stars is 
threaded with a magnetic field even though the hot B8V 
mass gaining star is not known to have a magnetic field. 
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